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基于 慰 夜 水 位 波动 法 估算 地 下 水 杂 散 发 量 的 研究 
一 一 以 河西 走廊 典型 绿洲 为 例 
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摘 要 : 使 用 多 种 昼夜 水 位 波动 法 (White 法 .Hays 法 .Loheide 法 ), 计 算 了 黑河 中 游荡 漠 绿 洲 过 渡 带 地 下 水 浅 埋 区 


(ET) AD. Wl tek HY AK ZSA (ET EAT ACEI To RRHH :在 几 种 算法 中 ,Hays 法 精度 最 高 ,其 次 是 White 法 ， 


生长 季 典 型 时 段 地 下 水 蒸 散发 (ET.) ,并 将 佑 算 结果 同 彭 曼 方法 获得 的 潜在 蒸 散 发 (PET) E-601 测量 的 水 面 蒸发 


Loheide 法 计算 精度 最 低 。 因 此 ,在 计算 逐日 ET: 时 可 优先 使 用 Hays 方 法 ,并 推荐 使 用 BE7" 来 检验 计算 精度 。 用 Lo- 


heide 法 计算 57 可 获得 较 高 的 精度 


(R=0.821 , P<0.01 ) ,但 具有 明显 时 滞 效 应 ,滞后 时 间 约 为 3 h。 这 些 计算 成 果 对 


当地 水 资源 的 合理 配置 与 可 持续 开发 利用 具有 一 定 的 参考 意义 和 应 用 价值 。 
关键 词 : 荒漠 绿洲 ; 地 下 水 ; 昼夜 水 位 波动 法 ; 地 下 水 补给 


绿洲 是 干旱 区 的 核心 地 理 单元 ,地 下 水 浅 埋 区 
(地 下 水 埋 深 <5~10 m) 受 地 下 水 变化 影响 显著 , 浓 
缩 了 地 下 水 灌溉 .土壤 水 盐 平衡 .生态 水 位 等 科学 
问题 ,是 绿洲 健康 指示 器 '"”。 地 下 水 是 地 下 水 浅 埋 
环境 中 植被 的 重要 水 分 来 源 , 植 被 提升 地 下 水 进行 
蒸 散发 (ET ) 是 干旱 区 地 下 水 重要 的 排泄 方式 。 
准确 估算 E7. 对 干旱 区 GSPAC( 地 下 水 -土壤 -植被 - 
大 气 -连续 体 ) 研 究 、 生 态 需 水 量 估算 .地 下 水 资源 
评价 具有 重要 意义 ,人 研究 结果 可 为 合理 利用 地 下 
水 、 遏 制 绿洲 生 态 环 境 恶 化 提供 科技 支撑 。 

影响 有 7 的 环境 因子 较 多 , 且 过 程 较为 复杂 , 实 
际 情况 下 很 难 对 其 进行 精确 估算", 现 有 的 ET. 估算 
方法 大 致 有 三 种 :(1) 水 量 平衡 法 :基于 水 量 平衡 原 
理 计算 7., 计 算 精度 受到 除 E7T, 以 外 的 其 他 水 均衡 
项 制约 ,日 影响 因素 较 多 ,很 难 实现 较 高 时 间 分 辨 
KWETE; (2) 极限 埋 深 法 :又 称 经 验 公 式 法 ,其 中 
经 验 系数 在 大 区 域内 应 用 时 误差 较 大 ,而且 此 方法 
受 环境 因子 制约 较 多 ,如 植被 盖 度 、 类 型 和 土壤 异 
质 性 等 5 ; (3) 昼夜 水 位 波动 法 :基于 地 下 水 水 位 罚 
夜 升降 波动 规律 来 计算 ET., 只 需要 给 水 度 和 地 下 
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水 位 数据 ( 埋 深 ) ,无 需 考 虑 区 域 土壤 异 质 性 和 地 
表 作 物 类 型 等 难以 精确 量化 的 信息 ,具有 所 需 参 
数 少 、 成 本 低 、 数 据 易 获 得 等 优势 ,所 以 在 近 些 年 
来 发 展 较 快 ” 

层 夜 水 位 波动 法 是 由 White 等 "于 1932 年 提出 
( 即 White 方法 ) ,此 方法 需要 满足 四 个 前 提 条 件 才 
可 行 :(1) 地 下 水 位 埋 深 的 下 降 是 由 植被 蒸 散发 吸 
收 地 下 水 导致 ;(2) 在 0:00 一 4:00 的 有 7 非常 微弱 ， 
EN ET,=0; (3) 0:00 一 4:00 的 地 下 水 补给 率 可 以 代表 
当日 所 有 时 间 的 补给 率 ;(4) 精确 计算 给 水 度 的 
值 。 有 很 多 学 者 针对 “四 个 假设 ”进行 了 改进 ,如 
Gribovszki 等 对 地 下 水 补给 率 和 给 水 度 的 改进 , 通 
过 水 力学 的 方法 计算 补给 率 , 使 其 计算 结果 更 加 接 
近 真 实 值 ;Meyboom'” 提出 使 用 “速效 给 水 度 ”(read- 
ily available specific yield ) E7440 24 7K BE ,并 推荐 给 
水 度 值 的 一 半 为 速效 给 水 度 ; Hays "针对 补给 率 的 
改进 ,将 水 位 变化 过 程 分 成 了 上 升 和 下 降 两 个 不 同 
时 段 ,分 别 计算 其 补给 率 与 消耗 率 , 进 而 估算 ET, 
Loheide 王 针对 给 水 度 和 补给 率 的 改进 ,引入 了 地 下 
水 位 县 加 原理 和 去 趋势 方法 ,并 给 出 了 给 水 度 计算 
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公式 ,可 以 精确 计算 小 时 分 辩 率 的 A7. ; Soylu 使 
用 傅 里 叶 变换 的 方法 对 地 下 水 位 数据 进行 拟 合 , 建 
立 87. 与 水 位 变 幅 间 的 定量 关系 ,该 方法 可 免 于 计 


于 南部 ,为 尾 殖 30 a 左右 的 新 绿洲 农田 及 种 植 年 限 
20 a 左右 的 防风 固沙 林 ( 图 1)。 该 地 区 属 典 型 的 半 
干旱 荒 并 草原 气候 ,年 平均 气温 7.6 %C ,降水 稀少 且 


算 补给 率 ,但 计算 过 程 更 为 复杂 。 目 前 ,昼夜 水 位 
波动 法 在 国内 干旱 区 应 用 较 少 ,研究 工作 尚 不 够 
深入 。 

本 研究 在 黑河 中 游 荒漠 绿洲 过 渡 带 地 下 水 依 
束 型 植物 覆盖 地 区 (Phreatophytes) 应 用 White , Lo- 
heide 和 Hays 等 方法 计算 ETs。 并 将 佑 算 结果 同 
PET、E7。 和 ET, 进行 相关 分 析 , 以 此 评价 各 方法 的 精 
度 及 适用 性 。 研 究 首 次 将 昼夜 波动 法 应 用 于 河西 
走廊 地 区 ,使 用 高 分 辨 率 数 据 ( 小 时 尺度 ) 计 算 87.， 
并 在 荒漠 绿洲 地 区 的 4 种 典型 地 表 景 观 进行 比较 分 
析 ( 农 田 `. 老 绿洲 、 防 护林 和 藉 漠 边缘 ) ,研究 成 果 可 
用 于 “一 带 一 路 ”区域 水 资源 的 评价 .荒漠 绿洲 地 区 
用 水 策略 的 改进 和 防治 土壤 盐 渍 化 提供 科学 依据 。 


1 研究 区 概况 


研究 区 位 于 甘肃 西部 .河西 走廊 中 段 的 黑河 流 
域 中 游 的 临 泽 县 车 漠 - 绿 洲 过 渡 带 (99"51'~100"30' 下 ， 
38°57'~39°42'N ,海拔 1304~2196 m) ,在 研究 区 内 农 
田 取 工 样 点 , 锁 龙 潭 湿地 南北 、 东 部 分 别 选 择 3 个 
样 点 :P 点 为 典型 老 绿洲 农田 ;Q 点 位 于 东部 , 属 其 
漠 边 缘 地 区 ,地 表 有 人 工种 植 的 权 梭 固沙 林 ;T 点 位 


图 1 研究 区 地 理 位 置 及 监测 


Fig. 1 Study area and monitoring well schematics 


季节 性 分 明 ,年 平均 降水 量 仅 110 mm ,年 潜在 蒸发 
量 >2000 mm, 干 时 指数 约 20, 生 态 环境 极其 脆弱 。 
地 势 南北 高 但 中 间 低 ,区 域内 地 下 水 主要 受 河道 补 
给 ,表层 土壤 水 分 主要 来 自 农 田 洪波 侧 向 渗 漏 和 浅 
层 地 下 水 补给 ”。 
研究 区 内 分 布 着 大 片 灰 棕 漠 土 ,境内 农作物 以 
玉米 为 主 , 降 水 稀少 ,农田 灌溉 基本 为 河道 水 与 地 
下 水 ,地 下 水 埋 深 较 浅 ( 埋 深 一 般 约 为 1~5 m) ,土壤 
盐 碱 化 现象 较为 严重 。 境 内 植被 生长 与 地 下 水 
关系 密切 ,多 为 地 下 水 依赖 型 植被 。 植 被 类 型 的 差 
异性 主要 是 由 地 下 水 的 埋 深 高 低 造 成 ,研究 区 内 主 
要 有 白杨 (Populus gansuensis) WX (Elaeagnus an- 
gustifolia) F2 (Caragana microphylla) 、 梭 梭 (Hal- 


oxylon ammodendron ) FU YY WI (Salix mongolica ) 等 时 


生 植物 ,与 新 垦 农 田 呈 镶嵌 分 布 格局 ws 。 
2 数据 与 方法 


2.1 数据 来 源 

(1) 水 位 埋 深 数据 是 在 研究 区 使 用 地 下 水 监测 
井 (PVC 管 做 的 简易 井 ,直径 约 $ cm) AY) HOBO JE 
力 传 感 器 来 获得 ,以 小 时 为 步 长 ,1h 测 1 次 。 监 测 


示意 图 
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项 目 包 括 大 气压 `. 地 下 水 埋 深 (mm) 和 水 温 , 对 水 位 
埋 深 数据 做 大 气 补偿 ,进而 得 到 真实 地 下 水 位 值 。 

(2) 气象 数据 采用 中 国 科 学 院 临 泽 站 标准 气象 
观测 场 的 实时 自动 监测 数据 ,时 间 分 辨 率 为 hn。 站 点 
位 于 甘肃 省 张掖 市 临 泽 县 五 里 墩 村 (100°07’42.0”E， 
39°20'59.0"N ) ,场地 由 固定 沙丘 平整 而 成 ,周围 植被 
以 梭 梭 、 树 柳 和 沙 束 为 主 。E-601 蒸 发 由 和 直径 20 
cm 的 Ga 蒸发 血 测 量 水 面 蒸发 量 , 时 间 分 辩 率 为 d。 

(3) 土壤 水 分 特征 曲线 :在 样 点 挖掘 土壤 剖面 ， 
使 用 标准 的 环 刀 (高 5 cm, 体积 110.11 em’) A 
埋 深 1.6 m 处 取样 品 ( 年 地 下 水 位 约 在 1.6m 人 处 上 下 
波动 ) ,然后 在 实验 室 中 使 用 H-1400PF 落 地 式 土壤 
用 高 速 冷冻 离心 机 测量 样品 的 土壤 水 分 特征 曲线 。 
2.2 ET, 估算 方法 
2.2.1 White 方法 (图 2) 

ET, =S,(24r +s) (1) 

SU ET, Ai Bk A AS a E (mm ) ;$S, 为 给 水 度 ;r 
为 0:00 一 4:00 AY HEF 7k Fb 29 (mm: he!) ;s AH F 
水 埋 深 日 变化 高 度 值 (mm) o 
2.2.2 Hays 方 法 (图 3) 两 个 假设 条 件 :(1) TAY Be 
ET.AL IS , BNET, =0; (2) 到 时 段 地 下 水 补给 率 可 
以 代表 当日 所 有 时 间 补 给 率 "”。 


H,-L 
ET, =S[(H,-L)+—— 
y T; 


式 中 :到 为 相 邻 两 日 地 下 水 埋 深 变化 过 程 线 的 第 1 
个 峰值 (mm); 厂 为 第 2 个 峰值 (mm) ;为 谷 值 ;7 为 
埋 深 下 降 时 段 ;为 埋 深 上 升 时 段 。 

2.2.3 Loheide 方 法 三 个 假设 :(1) 地 下 水 位 波动 
是 由 植被 吸收 地 下 水 进行 蒸 散 发 和 地 下 水 的 侧 向 


T] (2) 


图 2 White 方法 示意 图 
Fig.2 Diagram of White method 
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以 河西 走廊 典型 绿洲 为 例 


-850 — 地 下 水 埋 深 


图 3 Hays 方 法 示意 图 
Fig. 3 Diagram of Hays method 


补给 导致 ,不 考虑 其 他 因素 ;(2) FE 0: 00—6: 00 这 

WEN ET. AE FS SS , BILET,=0; (3) 区 域 地 下 水 埋 

深 的 变化 趋势 和 监测 井 埋 深 变化 趋势 一 致 "1。 
ET,=r(i)- st (3) 


式 中 :ET 的 量 纲 为 mm-h';r(1) 为 1 时 的 补给 率 ;WT 
为 地 下 水 埋 深 (mm) ;i 为 时 间 (h),dW7/di 为 埋 深 变 
化 率 (mm:h')。 
对 当日 至 翌日 0:00 一 6:00 地 下 水 埋 深 数据 做 
去 趋势 处 理 ,公式 为 ; 
WT = WT (krt + by) (4) 
式 中 ;WT 为 去 趋势 地 下 水 埋 深 (mm) ;如 为 去 趋势 
处 理 的 斜率 ;6 为 去 趋势 处 理 的 截 距 。 利 用 0:00 一 
6:00 的 WT, 得 到 相应 的 埋 深 变 化 率 ( wee ) , 建 
立 二 者 的 函数 关系 ， 
rwr,) e (5) 


RIET (Wor) KIHE H HABEAS E] Woe EV 

补给 率 , 公 式 为 : 
rQ)=S, [PWT r) + kl (6) 

继而 根据 公式 (3) 可 以 计算 出 7.。 
2.3 给 水 度 ($,) 估 算 

给 水 度 反 映 了 含水 介质 释 水 能 力 强 弱 。 据 Lo- 
heide EU E 2005 年 研究 发 现 ,当地 下 水 埋 深 <1.0 m 
时 ,Ss, 随 着 地 下 水 埋 深 而 变化 (水 位 深度 依赖 性 )。 
当地 下 水 埋 深 >1.0 m 时 ,S, 只 受 含水 介质 粒 径 结构 
(质地 ) 影 响 , 即 5, 会 越 来 越 接 近 某 一 个 常数 值 。 本 
人 研究 根据 土壤 水 分 特征 曲线 来 计算 土壤 水 分 特征 
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值 , 既 能 反映 出 地 区 土壤 特性 ,具有 一 定 的 代表 性 
又 避免 了 抽水 试验 计算 5, 的 复杂 不 确定 性 。 本 文 研 
究 区 地 下 水 埋 深 一 般 >1.0 m( 除 个 别 时 段 受 灌溉 影 
响 使 得 埋 深 <1.0 m, 本 文 不 予 考 虑 ), 所 以 本 文 认为 
S, 是 一 个 常数 值 ” 。 

可 根据 土壤 水 分 特征 曲线 进行 土壤 体积 含水 
量 (volumetric water content, VWC) 与 土壤 水 势 之 间 
的 换算 :土壤 水 势 为 0 时 VWC 为 饱和 含水 量 0., 土壤 
水 势 为 -1/3 bar 时 的 VWC 认 为 是 田间 持 水 量 0., 土 
壤 水 势 为 -15 bar 的 VWC 认 为 是 凋 萎 系数 9.” 2 。 

土壤 给 水 度 可 认为 是 饱和 含水 量 与 凋 萎 系数 
的 差 值 (土壤 最 大 有 效 水 量 )。 即 :5S,=0.-0.。 根 据 
Meyboom!"”! HY WEE ,本 文 计 算 时 使 用 速效 给 水 度 的 
值 (1/2 S,)。 
2.4.1 潜在 蒸 散 发 量 (PET) PET 是 指 在 植被 均匀 
履 盖 区 域内 ,地表 始终 属于 充分 供水 情况 下 的 蒸 散 
发 值 , 也 可 以 理解 为 地 下 水 埋 深 几乎 为 0 时 的 地 表 
蒸 散 发 。 本 研究 使 用 FAO 提出 的 改进 Penman- 
Monteith 计算 方法 。 


PET= 0408AR, -C)+ yg- gye lM.) e] (7) 
5 Aryl +0.34u,) 


和 温度 关系 曲线 的 斜率 (无 量 纲 ) ;及 ,为 地 表 净 辐射 
(MJ-m*-h");G Ay EA et (MJM h); ye 
RE Th BE BX (kPa + Co") ;7 为 单位 时 间 的 平均 气温 
(C) ;ww 为 小 时 时 段 内 的 平均 风速 (ms );e, 为 气温 


| 一 地 下 水 埋 深 


05-01 06-01 07-01 


为 四 时 的 饱和 水 汽 压 (kPa) ;e; 为 小 时 时 段 的 实际 水 
VUE (kPa) 5e.( The) -e WAKAAE (kPa), 
2.4.2 相关 性 分 析 Pearson 相关 性 分 析 方 法 ,作为 
统计 学 中 较为 常用 的 相关 性 分 析 方 法 ,Pearson 相关 
系数 能 够 有 效 简 洁 地 反映 出 两 种 不 同 要 素 之 间 相 
关 程 度 和 相关 方向 ,反映 出 两 要 素 之 间 的 相关 关系 
是 否 密切 中 ,用 R 表 示 , 计 算 公 式 如 下 : 


>o-ac 


人 Se. -7 


式 中 :R 为 相关 系数 (Pearson correlation coefficient) ; 
& 和 分 别 为 两 要 素 多 日 平均 值 ;a; 和 ;分别 为 要 素 
日 平均 值 。-1<R<1,R 的 绝对 值 越 大 表明 两 要 素 相 
关 性 越 强 ,R>0 表 明 两 要 素 呈 正 相 关 ,R<0 则 表明 两 
要 素 呈 负 相关 。 

在 研究 区 (地 下 水 浅 埋 区 ) 地 下 水 位 变化 与 地 
PARMA MIRA BA HY OY, PET 和 7 具有 颇 为 相似 
的 变化 波动 趋势 ,所 以 可 通过 对 PET 和 3 种 方法 计 
算 的 87. 之 间 进 行 Pearson 相关 分 析 , 并 做 显著 性 检 
验 , 用 来 检验 3 种 方法 的 计算 精度 。 


3 地 下 水 蒸 散发 量 计算 


31 地 下 水 动态 变化 

由 网 4 可 知 ,灌溉 发 生 后 地 下 水 位 迅速 上 升 , 而 
灌溉 停止 后 水 位 开始 下 降 , 但 下 降 过 程 较 为 缓慢 。 
在 研究 区 经 历 较 强 降水 事件 后 (2017 年 7 月 21 日 ， 


R= 


i=1,2,.…,n (8) 


0 
45 
+10 
T 
J 15 E 
420 名 
g 
J25 
| 30 
35 
40 
08-01 09-01 10-01 11-01 
日 期 /日 -月 


图 4 工 点 生长 季 降 水 和 地 下 水 埋 深 变化 过 程 线 


Fig. 4 Variation of precipitation and groundwater depth in L-point growing season 
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降水 量 约 19.2 mm) ,地 下 水 位 仍 未 产生 明显 响应 ， 
表明 降水 对 地 下 水 位 波动 的 影响 十 分 有 限 。 可 见 
在 研究 区 内 地 下 水 动态 变化 主要 受 灌 溉 等 人 为 因 
素 影响 ,而 降水 尤其 是 小 降水 事件 对 地 下 水 变化 的 
影响 并 不 明显 。 

考虑 到 昼夜 水 位 波动 法 的 适用 性 ,本 文选 择 时 
间 为 :2017 年 8 月 10 一 18 日 。 时 段 内 地 下 水 位 无 突 
变现 象 , 即 无 灌溉 .抽水 事件 发 生 。 时 段 内 仅 有 一 
次 3.6 mm 降水 事件 发 生 (2017 年 8 月 18 日 ), 且 降水 
有 件 发 生 在 16:00 以后。 由 于 小 降水 事件 产生 地 表 
径流 的 概率 极 低 , 本 文 忽 略 其 对 ET. 计 算 结 果 的 
影响 。 

根据 土壤 水 分 特征 曲线 计算 出 土壤 速效 给 水 
度 为 0.043( 表 1)。 

时 段 内 地 下 水 位 变化 昼夜 波动 显著 (图 5$) , 且 
整体 呈 下 降 趋 势 。 一 天 内 地 下 水 位 最 大 值 通常 出 
现在 8:00 一 10:00 之 间 ,最 小 值 通常 发 生 在 17:00 一 
19:00 之 间 。 

3.2 ADIZ ET, 

受 计算 方法 的 局 限 性 ,本文 计 算 ETKA 
9 d。 为 了 确保 Pearson 相关 分 析 的 样本 数量 充足 且 
具有 代表 性 ,本 研究 使 用 Kolmogorov-Smirnov 法 (K- 
S 检 验 ) 检 验 了 样本 数据 的 正 态 性 ,结果 表明 3 种 方 


alin 


R1 工 点 土壤 水 分 特征 值 及 给 水 度 


Tab. 1 Characteristic value of soil moisture and specific 


yield at L point 
特征 值 ”饱和 含水 量 田间 持 水 量 调 蓉 湿度 给 水 度 速效 给 水 度 
土壤 水 势 0 U3bar 15 bar 
VWC 0.352 0.346 0.267 0.085 0.043 
-800 
-900 
g -1000 
Rè -1100 
| -1200 
5 -1300 
_1400| 一 地 下 水 埋 深 
DNS 闪光 
BLK KKK KK HH 
日 期 /月 -日 
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以 河西 走廊 典型 绿洲 为 例 


法 计算 的 有 7, 均 服从 正 态 分 布 (P<0.05)25 。 

利用 White 和 Hays 法 计算 获得 的 该 时 段 内 的 
ET.AAA 6~9 mm- d'( 16), BET, 49 PET ET) All ET, 
变化 趋势 一 致 ,但 由 Loheide 法 计算 的 Z7. 在 8 月 11 
日 和 8 月 17 日 出 现 异常 ,甚至 出 现 了 ET.(11.2 mm: 
d!)>PET(7.7 mm: dO HITED 

HX 2 aA, AS EA ERE H ET, 5 PET, 


12 


AAR (mm: d’) 
a oo 


—-®—Loheide 
== ET, 


一 = 一 White 
2 —s— PET 


~-a--- Hays 
----- ET, 


0 
8-10 8-11 8-12 8-13 8-14 8-15 8-16 8-17 8-18 
日 期 /月 -日 
图 6 LARERE 


Fig.6 L-pointevapotranspiration change chart 


ET), ET, Z lil K 4H OG YE : Hays>White>Loheide. 
此 ,Hays 方 法 在 计算 日 分 辩 率 ET, BR Be HF (R> 
0.8) ,White 法 其 次 ,Loheide 法 最 差 。 验 证 各 方法 计 
算 精度 最 好 使 用 E-601 测 得 的 £7, 数据 ,其 次 是 Bw 
测 得 的 £7 数据 和 PET。 
3.3 小 时 分 辨 率 ET， 

由 于 Loheide 法 计算 的 87. 与 PET 相 关 性 较 弱 ， 
为 剖析 其 原因 ,本 研究 使 用 Loheide 方 法 计算 了 小 时 
AY PER ET,, 并 与 逐 时 PET 进 行 相关 分 析 。 结 果 表 
HH :其 相关 性 仍然 较 弱 但 依旧 显著 (R=0.463, P< 
0.01) 。 通 过 对 PET 时 滞 处 理 后 ,相关 性 明显 提高 
( 表 3) , 且 将 时 滞 设 置 为 2~4h, 相 关 性 较 好 (有 
0.75, P<0.01) , 当 设 置 时 沾 为 3h 时 ,相关 性 最 好 (R= 


表 2 ET,5 PET. ET AN ET, A Pearson 相关 分 析 


Tab. 2 Pearson correlation analysis of ET, with 


Fig. 5 Groundwater burial depth change chart in L-point 


图 $ 工 点 研究 时 段 地 下 水 埋 深 变化 


PET, ET, and ET, 
PET ETo ET, 
White 0.761" 0.792" 0.739° 
Loheide -0.105 0.005 0.001 
Hays 0.843” 0.897” 0.872" 
注 :** 在 0.01 级 别 ( 双 尾 ) ,相关 性 非常 显著 ;* 在 0.05 级 别 ( 双 尾 )， 


相关 性 显著 。 下 同 。 
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表 3 各 样 点 时 ET 与 时 PET 的 Pearson 相关 分 析 
Tab.3 Analysis results of Pearson correlation between 


hourly ET, and hourly PET at different sample 


趋势 地 下 水 位 埋 深 (7 与 PET,ET 这 3 种 信号 均 
具有 明显 的 极 值 ,但 W7T 信 号 的 峰值 区 域 更 加 圆 
整 。 这 3 个 信号 峰值 和 谷 值 在 时 间 上 并 没有 完全 重 


县 , 出 现 极 值 时 间 前 后 :Tow>E7.>PET,7Tw 最 先 达 
到 峰值 ,其 中 PET 信 和 号 极 值 滞后 7 约 6h,PET 信 
号 极 值 济 后 11. 约 3 h, 佐 证 了 上 述 滞后 时 长 为 3h 
左右 的 观点 。 

图 7 和 图 8 中 有 部 分 时 间 计 算 的 E7. 是 负 值 ,这 
可 能 是 因为 压力 传感器 (测量 地 下 水 埋 深 ) 产 生 的 噪 
声 使 得 去 趋势 水 位 埋 深 变化 率 dW7wwdi 5T (WTor)2 


样 点 È P Q T 
时 间 08-10— 06-24— 07-06 一 07-10 一 
08-18 06-30 07-11 07-17 

无 时 滞 0.463" 0.380" 0.549" 0.587" 
Ii lh 0.649" 0.585" 0.690” 0.738” 
Hih 0.772” 0.722" 0.784" 0.825” 
Wi 3h 0.821” 0.783“ 0.822” 0.833” 
INir4 h ”0.795™ 0.785" 0.771" 0.757" F 
rSh ”0.702™ 0.714" 0.650” 0.616” 
时 滞 6h ”0.563” 0.593" 0.459” 0.437” 


È IRR AER HFT RT SE AB 


0.821, P<0.01). 

对 各 样 点 计算 的 E7. 与 PET 数 据 进 行 了 对 比 ， 
结果 表明 在 不 同 植被 覆盖 条 件 下 均 存在 不 同 程度 
时 清关 系 。Q 两 点 时 滞 为 3 hp,P 点 时 涉 为 4h( 表 
3) ,时 滞 处 理 后 的 E7, 和 PET 相 关系 数 尺 明显 提高 ， 
RMSE 明显 降低 ,使 用 Loheide 法 计算 时 ET 效果 较 
好 (图 7)。 可 见 , 在 黑河 流域 中 游 地 区 通过 地 下 水 
位 波动 法 获得 的 ET, UE aed FE LE I ie PR, Air Jia 
间 在 3 左右 。 

Loheide 法 能 够 较 好 地 利用 高 频 地 下 水 位 数据 
(分 辨 率 为 1h 或 者 更 高 ) 来 计算 时 7,; 与 逐 时 ET, 
相 比 ,PET 值 通常 会 延迟 几 个 小 时 (时 滞 性 )。 

地 下 水 埋 深 越 浅 表明 对 地 表 植 被 供水 越 充 分 ， 
说 明 有 7 与 地 下 水 埋 深 关 系 密切 。 从 图 8 可 看 出 ,去 


种 信号 不 匹配 导致 的 结果 '"。 


4 计算 成 果 合 理性 分 析 


昼夜 水 位 波动 法 在 黑河 流域 中 游 地 区 具有 很 
好 的 适用 性 ,估算 的 1, 值 较 大 (6~9 mm.d) ,这 一 
结果 与 孙 海 涛 等 ”在 中 国 北方 其 他 干旱 AB 
地 区 (毛乌素 沙漠 ,鄂尔多斯 ) 的 研究 结果 相 接 近 ， 
结果 可 信 。 

4.1 验证 计算 精度 数据 选择 

由 表 2 可 知 , 验 证 3 种 方法 的 计算 精度 最 好 使 
用 ET 数据 (日 分 辨 率 E7.)。 这 可 能 与 仪器 自身 结 
构 有 关 ,E-601 的 直径 较 大 ,安装 结构 更 加 合理 科 
学 ,其 所 测 得 的 水 面 蒸 发 量 E7, 更 具 代 表 性 。 而 Bw 
蒸发 耿直 径 只 有 20 cm, 并 不 能 很 好 的 代表 区 域 水 
面 蒸 发 能 力 。PET 是 依据 芯 曼 公式 计算 出 的 结 
属于 模拟 数据 并 非 实 测 数据 , 且 彭 曼 公 式 属于 半 经 


2; R=0.463 è LÄ 2; R=0.380 © PA R=0.549 2j R=0.587 
z RMSE-0343 拟 合 线 Een RMSE=0.295 一 一 拟 合 线 = 2} e RMSE=0.341 p RMSE=0,318 
T ® oos 7 . 了 oe 4 s I" au 
= 1 < 1 z a vay 
: g 中 二 Tite x 
Š = = E | Peat 
>p p 0 = RY 0 a" 
四 z a OF ig n" T 
a 。 d 4 ”一 拟 合 线 
00 05 10 1.5 00 05 1.0 0.0 0.5 10 15 0.0 0.5 10 15 
PET/(m-h”) PET/(m:h”) PET/(m:h”) PET/(m-h”) 
2f R=0.821 27 R=0.785 R=0.822 2f R-0833 ， 
= | RMSE=0.203% © | RMSE=0.197 „e > 2} RMSE=0.233° 个 “| RMSE-0.218 ， 
| ee ee x] ° =n 
. 1 C] . rn . “ee p 1 
g oe ge E [ 
E i E El yarn E 
图 0 S So Es 8 0| apet 
Fats 。 LETH) wota Qua (RY Fir) oe TAITH) 
E š 拟 合 线 si 一 一 拟 合 线 e 一 一 拟 合 线 
00 05 10 1.5 00 05 1.0 0.0 05 10 15 0.0 05 10 15 
PET/(m-h*) PET/(m-h*) PET/(m-h*) PET/(m:h*) 


图 7 WLU ke HS AR Ja AY ET, A PET 的 散 点 关系 
Fig. 7 Contrastive scatter plots of four-sample time-lag ET, and PET 
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蒸 散发 (mm'h-) 
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WTpr/mm 


| —80 


| | | 
08-10 08-11 08-12 08-13 08-14 
日 期 /月 -日 


| L 1 | 
08-15 08-16 08-17 08-18 08-19 


图 8 WTor, BT. 与 PET 随 时 间 变 化 关系 


Fig. 8 WTor, ET,and PET time-varying chart 


验 半 理论 公式 ,理论 依据 不 够 充分 等 原因 ,所 以 ,将 
PET 作 为 验证 数据 也 具有 一 定 的 局 限 性 。 
4.2 ET. 与 PET 的 滞后 关系 

Loheide 法 能 够 满足 计算 小 时 分 辨 率 ET, Ee 
计算 更 高 分 辨 率 E7., 时 清 处 理 后 的 计算 精度 明显 
提高 ,RMSE 值 也 明显 降低 。 乔 将 Loheide 法 应 用 于 
其 他 类 似 地 区 ,要 考虑 时 渍 效应 。 混 后 时 间 长 短 可 
能 与 当地 水 文 、 气 象 .土壤 质地 和 植被 类 型 等 因素 
相关 ,在 本 研究 区 不 同类 型 生态 系统 存在 的 清 后 时 
长 约 为 3 h( 地 表 履 盖 梭 梭 林 的 老 绿洲 生态 系统 滞后 
时 长 约 为 4h)。 

时 涡 现 象 可 能 是 由 于 植物 的 水 分 运输 机 制造 
成 的 :因为 植物 吸水 耗 水 生理 过 程 ,地表 蒸 散 发 与 
地 下 水 埋 深 的 响应 过 程 不 可 能 同时 进行 ,滞后 的 时 
间 可 能 是 植被 提前 从 地 下 水 中 吸收 水 分 运 移 到 冠 
层 进行 储存 ,而 后 进行 蒸 散 发 这 一 过 程 所 需 时 间 。 
地 表 植 被 冠 层 也 可 储存 大 量 水 分 ,并 且 这 种 水 储量 
也 能 够 反映 出 植被 燕 散 发 与 根系 吸收 地 下 水 或 土 
二 水 之 间 存 在 的 时 间 涡 后 现象 下。 此 外 , 据 L[i 等 " 
研究 结果 也 表明 ,土壤 水 分 在 非 饱 和 带 运 移 ( 上 升 ) 
过 程 中 通常 也 会 有 3 ~ 4 h 的 滞后 现象 。 

4.3 不 确定 性 分 析 

(1) 计算 方法 的 不 确定 性 :Hays 法 是 根据 地 下 
水 埋 深 变化 过 程 线 的 峰 谷 值 选 择 补 给 阶段 及 耗 散 
阶段 ,避免 了 人 为 选择 ET, =0 的 时 段 导 致 计算 的 不 
确定 性 ,更 接近 于 真实 情况 ,所 以 此 方法 计算 精度 


Rei), White ITY ET, PET Fil ETE fb 
相关 ,R>0.7(P<0.05)。 但 由 于 人 为 选择 ET,=0 的 时 
间 段 具有 一 定 的 局 限 性 ,其 计算 精度 略 低 于 Hays 
法 ,表明 White 法 是 一 个 较为 经 典 的 计算 方法 。 而 
Loheide 法 计算 的 87. 与 PET、ET, 和 ET 几乎 毫 无 关 
系 。 这 可 能 与 Loheide 法 自身 缺陷 有 关 , 采 用 去 趋势 
手段 削弱 了 原始 地 下 水 埋 深 的 部 分 信号 值 ,掩盖 了 
一 些 有 用 信息 。 此 外 假设 £7.=0 的 时 段 是 0:00 一 
6:00, 时 有 段 代 表 性 不 足 , 且 日 分 辩 率 数据 本 身 尺度 
较 大 ,对 计算 结果 也 有 一 定 影响 。(2) 布设 .检查 仪 
器 和 传 感 硕 自动 监测 数据 均 对 数据 的 采集 产生 影 
响 ( 噪 声 ) ,模拟 数据 (PET) 和 实测 数据 均 包 含 噪 声 ， 
对 计算 结果 有 影响 呈 。(3) 地 下 水 埋 深 的 大 小 以 及 
埋 深 的 变化 速率 直接 影响 土壤 5,, 埋 深 也 在 时 刻 变 
化 着 ,具有 明显 的 昼夜 波动 现象 。 当 水 位 埋 深 的 上 
升 或 下 降 速率 很 大 时 ,$, 并 不 是 恒定 的 ,但 在 本 文 的 
计算 过 程 中 认为 其 是 定 值 ,此 时 引入 了 瞬 态 误差 。 
(4) 选择 典型 的 研究 时 段 要 考虑 到 无 较 大 降水 事 
件 .无 水 和 泵 抽取 地 下 水 以 及 农田 灌溉 等 大 型 抽水 或 
补给 事件 发 生 ,这 些 因素 制约 了 计算 时 段 的 选择 ， 
使 得 符合 假设 条 件 的 原始 样本 数据 较 少 ,对 计算 结 
果 也 会 产生 一 定 的 影响 。 


5 结论 


通过 上 述 分 析 计算 ,得 到 以 下 结论 : 
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(1) White, Loheide 和 Hays 三 种 方法 在 河西 走 
廊 中 游 地 区 均 具 有 一 定 适用 性 ,其 中 Hays 法 是 根据 
地 下 水 埋 深 变化 过 程 线 的 特性 自然 选择 补给 和 耗 
散 时 段 ,更 接近 于 真实 情况 ,计算 效果 最 佳 ,而 Lo- 
heide 和 White 法 是 人 工 选 择 补给 和 耗 散 时 段 ,具有 
一 定 的 局 限 性 ,推荐 使 用 Hays 法 计算 日 尺度 ET,。 

(2) Loheide 法 可 以 计算 更 高 分 辨 率 E7,, 如 15 
min 分 辩 率 。 给 水 度 的 值 会 直接 影响 到 计算 结果 ， 
选择 合适 的 给 水 度 十 分 重要 。 

(3) 若 将 Loheide 法 应 用 于 其 他 相似 地 区 , 需 
虑 滞后 时 间 效 应 。 在 本 研究 区 PET 落 后 有 7 的 滞后 
时 长 约 为 3h。 

(4) 使 用 结构 合理 和 安装 科学 的 E-601 节 发 耿 
实测 的 57, 数 据 验证 三 种 方法 计算 日 7 的 精度 效 
果 最 好 。 
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Abstract: Daily groundwater evapotranspiration (ET,) was estimated for the shallow groundwater area of the 
new oasis in the middle reaches of the Heihe River basin using diurnal water level fluctuation methods (White, 
Hays, Loheide). The results were compared with the potential evapotranspiration (PET) estimated using Penman 
method, water surface evaporation (E7)) measured using E-601 pan, and water surface evaporation (£7') measured 
using Pv, and the differences between ET, and PET (ET, and ET) were analyzed using Pearson correlation 
coefficient. The results showed that among these calculation methods, Hays method performed the best with the 
highest accuracy, followed by the White method, and the Loheide method performed the worst with the lowest 
accuracy. Therefore, Hays method is recommended for estimation of daily ET, and ET) is recommended for 
testing the estimation accuracy. Loheide method showed high ET, accuracy (R=0.821, P<0.01) with approximate 
3 hr time lag, which was also found in other areas of the Heihe River basin (old oasis, desert edge, shelterbelt). 
The results obtained in this study will be helpful and may offer scientific support for rational allocation, 
sustainable development, and utilization of regional water resources in the Heihe River basin. 


Keywords: desert oasis; groundwater level; diurnal groundwater variation; groundwater recharge 


